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NOMENCLATURE 



A - constant in regression rate equation (7) 

B - blowing parameter 

Cp - specific heat at constant pressure 

E - activation energy 

h - enthalpy 

h° - heat of formation 

k - turbulence kinetic energy 

m - mass flow rate 

R - average molecular weight 

m^ - mass fraction of species i 

p,P - pressure 

- heat transfer flux 
r - radial distance 

r - fuel regression rate 

R - universal gas constant 

R.j - reaction rate of species i 

T - temperature 

u - axial velocity 

v - radial velocity 

w - tangential velocity 

r. - effective transport coefficient of species i = y/a. 

<5 - incremental distance from wall 

e - turbulence dissipation rate 

u - effective viscosity 

p - density 

a.. - Prandtl or Schmidt number for species i 

fraction (mass fraction of total fuel) 

vi i 



mixture 



subscri pts 



a 


- 


ai r 


c 


- 


chamber 


CO 


- 


carbon monoxide 


C0 2 


- 


carbon dioxide 


CH 


- 


intermediate hydrocarbon 


f 


- 


solid fuel 


fu 


- 


unburned fuel vapor 


g 


- 


within the fuel grain 


H 2 


- 


hydrogen 


h 2 o 


- 


water 


°2 


- 


oxygen 


P 


- 


near wall node point, poi 


T 

w 


- 


wall temperature 


rad 


- 


radiation 


w 


_ 


fuel wall 
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I. INTRODUCTION 



The solid-fuel ramjet (SFRJ) most often consists of a solid-fuel grain 
which provides the walls for the combustion chamber. A sudden expansion at 
the air inlet (either axial or side-dump) can be used to provide flame 
stabi lization. Downstream of the flow reattachment a turbulent boundary layer 
develops and includes a diffusion-controlled flame between the fuel-rich zone 
near the wall and the oxygen-rich central core. Due to that diffusion flame, 
heat is transferred by convection and radiation to the solid surface, causing 
vaporization of the fuel. 

At the Naval Postgraduate School both mathematical modeling* -1 * and 
experimental efforts^-lO have been conducted to determine the effects of 
design and operational variables on the obtainable performance. 

The computer model simulation of the SFRJ combustion process has evolved 
from an original stream-function vorticity formulation* to a 
primiti ve-variable (pressure, velocity) model which includes an aft mixing 
chamber. ^ These axisymmetric models have not included radiative heat transfer 
to the fuel surface. This has limited the utility of the model, since many 
al 1 -hydrocarbon fuels produce significant amounts of radiative transfer 
through the generation of carbon particulates in the flame zone. Recently*, a 
simplified radiative model was included in the primitive-variable model and 
resulted in improved prediction of the fuel regressi on-rate profile. The 
primitive-variable codes were based on the Champion 2/E/FIX computer program 
developed by Pun and Spalding**. The Champion code is a two-dimensional code 
and therefore, could not be used to model two important SFRJ geometries; 
side-dump and bypass. The side-dump configuration has air injected through 
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the side of the fuel grain and the dome region is constructed of fuel. The 
bypass configuration ducts a portion (typically 50%) of the total air flow 
into an aft mixing chamber. A full three-dimensional (3-D) model is needed to 
predict the flow-field in these configurations. 

The side-dump configuration is a very difficult one to develop/optimize 
using only experimental data. There is a complex interaction between the 
location/size of the inlet dump, the length of the dome region upstream of the 
dump, and the resulting fuel regression rate pattern. Regression rate data 
for the side dump geometry have not appeared in the open literature. However, 
several general observations can be made. Swirling the inlet flow can 
increase and smooth-out the downstream fuel regression rate pattern. 

Regression rates within the dome region, upstream of the dump, can increase or 
decrease, depending upon the dome length, number of inlets and amount of 
swirl. 3-D models are, in this case, a necessity in the combustor development 
process . 

Several 3-D, elliptic, primitive-variable codes have been developed. . One 
such code, developed by the AiResearch Manufacturing Company, 12, 13 -j s an 
extension of the Champion code. The former computer code was developed to 
model gas turbine combustors and formed the basis of the 3-D, SFRJ model 
discussed herein. 
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II. MODEL OVERVIEW 



The Garrett/AiResearch model 
and includes finite-rate chemical 
u» v, w, h, k, £, $» * m CH ’ n H 



is 3-D, elliptic, for steady, subsonic flow 
kinetics. The dependent variables are 
and In addition, if radiation is 



included (a six-flux model), then three radiation flux vectors are calculated. 
Soot emissions can be calculated if desired. Liquid spray calculations can 
also be made but were not necessary in the SFRJ model. The turbulence model 
is the same as used in the Champion program, namely the modified two-equation 
(k-e) model of Jones, Launder and Spalding^ - ^. 

Hydrocarbon combustion is considered to be a four-step process as 
f ol 1 ows : 



CH + C H 0 + H 0 
x y x y-2 2 

C H 9 + £ 0 o - xCU + ^ H, 
x y-2 22 22 



CO + ~ 0 2 + C0 2 



H 2 + J 0 2 H 2 0 



( 1 ) 

( 2 ) 

(3) 

(4) 



Arrhenius rate expressions are used for each of the above reactions. These 
chemical kinetics limited rate expressions were compared to rate expressions 
based on an eddy-breakup model. The latter expressions attempt to account for 
the effects of turbulence intensity and scale and species concentration on the 
consumption rate of a particular reactant. The two-part boundary layer 
(divided by y + = 11.5) used in the Champion code is also employed. Details of 
the computer program, the differential equations and the line by line, finite 
difference solution procedure are presented in references 12 and 13. 
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III. MODEL MODIFICATIONS AMD SOLUTION PROCEDURES 



The major modifications required of the Garrett/Ai Research model were 
geometric and the inclusion of a fuel wall on which finite-rate chemical 
kinetics could occur. The geometric changes required to model the SFRJ 
configurations were readily incorporated. Several additional variables and 
subroutines were necessary to include the regressing fuel surface effects. In 
addition, fuel property subroutines were changed to reflect the properties of 
Plexiglas (PMM) and HTPB, the two solid fuels used in the present study. 

The governing equations on the boundary^ for unity Lewis number are: 



Energy: 



q = (pv) (h -h.) + r. (|^) 
rad v 'w' w f ' h '3r ; w 

w 

(q r ^ d can be calculated within the existing model 
neglected in the present investigation.) 



( 5 ) 

but was 



3 m. 

Species: R. =T. (^p-) w + (pv) (m. ) (6) 

w w w g 

(the values of Rj were calculated in the same manner as for the 
internal nodes) 

It was assumed that the fuel regression rate could be represented in an 
Arrhenius format:^ 



P f r = (pv) w = A Pf e" E/RT w (7) 

Using the Couette flow approximation for the boundary layer behavior with mass 
transf er^ > ^ 



r 



*n(l+B) 
B=0 B 



( 8 ) 



where B can be evaluated from the enthalpy as 
h -h 

D = P w 

h -h, - (q "./(pv) )" 
w fg v ^rad ' ' w ' 



( 9 ) 
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The enthalpy can be expressed as 



l = E m ( (h” + J " C dT) 

w i' f i P, 

ref 



( 10 ) 



The solution procedure for the boundary conditions on the fuel wall was 
as follows: 

(a) estimate T and B 

(b) calculate (pv) w from (7) 

(c) calculate r. and T h from (8) 

w w 

(d) calculate n,. for all species from (6) 

w 

(e) calculate the other species (m u nJ m nv , m„ n , n.. ) from 

species conservation equations 2 22 

(f) calculate h w from (10) 

(g) calculate B from (9) with = 0 

Iterate (c) to (g) until B converges 

(h) calculate (pv) from (5) 

w 

(i) calculate T w from (7) 

Iterate (b) to (h) until (pv) w converges 

(j) calculate p w from 



w 



= P 
~T* 



(k) calculate v from ^ pV ^w 
v ' w 

p w 
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It was found, as previously found by Metzer* and Stevenson and Netzer^, 
that the near-wall turbulence dissipation had to be increased on the step face 

(e^=kp3/2/o.4<5 ) and that the grid spacing adjacent to the fuel surface had to 

be maintained slightly less than that required for y + =11.5 in order to 

P 

obtain temperature distributions in qualitative agrement with experimental 
data. 

In order to obtain convergence it was necessary to input initial 
conditions which were in qualitative agreement with the final solution (for 
example, the fuel mass fraction had to be specified as decreasing from the 
wall) and under relaxation (typically 0.1 to 0.5) had to be used. 

The accuracy of the solution for the entire combustor of the solid fuel 
ramjet is particularly sensitive to the boundary conditions on the fuel wall. 
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IV. RESULTS AND DISCUSSION 



Four geometries were investigated (Fig. 1); axi symmetri c with axial 
inlet, bypass with axial upstream inlet , single side-dump and dual side-dump. 
All cases included an aft nixing chamber. For the axi symmetri c , axial-inlet 
geometry the predicted reattachment "points" (where u near the wall was 
interpolated to zero) are shown in Fig. 2. The reattachment "points" were in 
good agreement with the maximum heat transfer locations measured by Krall and 
Sparrow^ but shorter than those measured by Phaneuf and Netzer® in 
nonreacting flow with wall mass addition. 

Normalized regression rate profiles are shown in Fig. 3. The 
experimentally measured profiles for PMM fuel were very sensitive to inlet 
turbulence/distortion. However, the 3-D code (with finite rate kinetics) 
predicted results very nearly identical to the mixing limited, 2-D code^, 
i.e., the maximum regression rate was predicted to occur upstream of the 
experimentally measured point. The average regression rate was in good 
agreement with experiment. 

Examples of typical radial profiles for u, m n , m M , and T 

2 2 

near the grain exit are shown in Fig. 4. The effects of including finite-rate 
kinetics and/or the injection of 03 that occurs into the boundary layer near 
the reattachmnt zone are readily apparent; a broad zone of high temperature 
near the maximum value and finite concentrations of oxygen below the flame. 

Figure 5 shows regression rate profiles at three axial locations for the 
single side-dump geometry without inlet air swirl. Non-uniform, high 
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regression rates occur in the dump region. Lower regression rates in the dome 
region and more uniform regression rates near the grain exit are also 
observed. 

Figure 6 shows the predicted regression rate results with two, 180° 
opposed side-dumps, with and without inlet air swirl. The effects of inlet 
swirl on this particular geometry are readily apparent. Inlet swirl decreased 
the regression rates in the dome, but significantly increased the rates 
downstream. Regression rate uniformity was also significantly improved except 
adjacent to the inlet dump. 

Figure 7 shows the results of increasing the dome length and/or the inlet 
dump area. Increasing the inlet area resulted in increased fuel consumption 
in the dome region and increased nonuniformity downstream. Lengthening the 
dome above the reference value significantly reduced head-end fuel regression 
rates without major effects in the downstream regions. 
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V. CONCLUSIONS 



The model predictions show that significant variations in fuel 
regression rate (both circumferentially and axially) are caused by the manner 
in which the air is injected into the fuel grain. In general, inlet swirl 
increases and smooths-out the fuel regression rate downstream of a side-dump 
inlet. The effects of inlet location and size on the regression rate behavior 
are difficult to estimate a-priori for a specific set of test conditions. The 
model must be used to predict the expected results. Much additional work is 
required in order to validate the predictions of the 3-D SFRJ model. However, 
the ability to examine the effects of inlet air location and flow 
characteristics on the fuel regression rate pattern (and the combustion 
efficiency, etc.) has provided a needed improvement to the earlier 2-D model. 

Current work is being directed at incorporation of the soot 
production/consumption and radiation subroutines into the SFRJ configuration 
and combustion environment. Regression rate profiles for the more recent 
configurations are also being compared to experimental data as they become 
avai lable. 
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Figure 1. Solid Fuel Ramjet Geometries in Model. 
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Figure 2. Reattachment Locations for Axial-Inlet Flow. 
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Figure 3. Fuel Regression Rate Profiles for Axial-Inlet Flow. 
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Figure 4. Predicted Radial Profiles for Gas- Properties Near the Grain 
Exit, Axial-Inlet Flow. 
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Figure 5. Fuel Regression-Rate Profiles, Single-Side Dump, No Inlet Swirl. 
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Figure 6. Effects of Inlet Swirl on Fuel Regression Rate Profiles, 180° Opposed 
Side- Dumps . 
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Figure 7. Effect of Inlet Area and Dome Length on Fuel Regression Rate Profiles, 180 
Opposed Side-Dumps with Inlet Swirl. 
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